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We have searched for Gravitational Waves (GWs) associated with the SGR 1806 — 20 hyperflare 
of 27 December 2004. This event, originating from a Galactic neutron star, displayed exceptional 
energetics. Recent investigations of the X-ray light curve's pulsating tail revealed the presence of 
Quasi-Periodic Oscillations (QPOs) in the 30 — 2000 Hz frequency range, most of which coincides 
with the bandwidth of the LIGO detectors. These QPOs, with well-characterized frequencies, can 
plausibly be attributed to seismic modes of the neutron star which could emit GWs. Our search 
targeted potential quasi-monochromatic GWs lasting for tens of seconds and emitted at the QPO 
frequencies. We have observed no candidate signals above a pre-determined threshold and our 
lowest upper limit was set by the 92.5 Hz QPO observed in the interval from 150 s to 260 s after 
the start of the flare. This bound corresponds to a (90% confidence) root-sum-squared amplitude 
h?ss-det = 4-5 x 10~ 22 strain Hz~ 1//2 on the GW waveform strength in the detectable polarization 
state reaching our Hanford (WA) 4 km detector. We illustrate the astrophysical significance of the 
result via an estimated characteristic energy in GW emission that we would expect to be able to 
detect. The above result corresponds to 7.7 x 10 46 erg (= 4.3 x 10 _8 Mqc 2 ), which is of the same 
order as the total (isotropic) energy emitted in the electromagnetic spectrum. This result provides 
a means to probe the energy reservoir of the source with the best upper limit on the GW waveform 
strength published and represents the first broadband asteroseismology measurement using a GW 
detector. 

PACS numbers: 04.80.Nn, 07.05.Kf, 95.85.Sz, 04.30.Db, 95.55.Ym, 04.40.Dg, 97.60.Jd, 97.10. Sj 



I. INTRODUCTION 

Soft Gamma-ray Repeaters (SGRs) are objects that 
emit short-duration X and gamma-ray bursts at irregu- 
lar intervals (see [l| for a review). These recurrent bursts 
generally have durations of the order of ~ 100 ms and 
luminosities in the 10 39 — 10 42 erg/s range. At times, 
though rarely, these sources emit giant flares lasting hun- 
dreds of seconds (see for example [USUI) withpeak elec- 
tromagnetic luminosities reaching 10 47 erg/s [5|. Pulsa- 
tions in the light curve tail reveal the neutron star spin 
period. 

Quasi-Periodic Oscillations (QPOs) @, i, B 03 in 
the pulsating tail of giant flares were first observed for the 
27 December 2004 event of SGR 1806-20 by the Rossi X- 
Ray Timing Explorer (RXTE) and Ramaty High Energy 
Solar Spectroscopic Imager (RHESSI) satellites [f| 0, H|. 
Prompted by these observations, the RXTE data from 
the SGR 1900 + 14 giant flare of 27 August 1998 was 
revisited [ll|. Transient QPOs were found in the light 
curve pulsating tail at similar frequencies to the SGR 
1806 — 20 event, suggesting that the same fundamental 
physical process is likely taking place. 

Several characteristics of SGRs can be explained by the 
magnetar model fl2l ]. in which the object is a neutron 
star with a high magnetic field (B ~ 10 15 G). In this 
model the giant flares are generated by the catastrophic 
rearrangement of the neutron star's crust and magnetic 
field, a starquake [l3l fl4|. 

It has been suggested that the star's seismic modes, 
excited by this catastrophic event, might drive the ob- 
served QPOs @, 0, B Il5j . which leads us to investigate 
a possible emission of Gravitational Waves (GWs) asso- 
ciated with them. There are several classes of non-radial 
neutron star seismic modes with characteristic frequen- 
cies in the ~ 10 — 2000 Hz range [pq . Toroidal modes of 
the neutron star crust are expected to be excited by large 
crustal fracturing (see [1, 0, H Qj}), though these modes 



may be poor GW emitters. However, crust modes could 
magnetically couple to the core's modes, possibly gen- 
erating a GW signal accessible with today's technology 
(see [IHHilHil). Other modes with expected frequencies 
in the observed range are crustal interface modes, crustal 
spheroidal modes, crust/core interface modes or perhaps 
p-modes, g-modes or f-modes. The latter should, in the- 
ory, be stronger GW emitters (see for example [2ll. I22I]). 

In addition, it has been noted [23[ that a normal neu- 
tron star can only store a crustal elastic energy of up to 
~ 10 44 erg before breaking. An alternative to the con- 
ventional neutron star model, that of a solid quark star, 
has also been proposed in several versions [23, [H, [25|, [2(| . 
In this case an energy of ~ 10 46 erg (as observed for this 
flare) is feasible, and thus the mechanical energy in the 
GW-emitting crust oscillations could be comparable to 
the energy released electromagnetically. This was also 
noted by Horvath [27|, who in addition estimated that 
LIGO might be able to detect a GW burst of compara- 
ble energy to the electromagnetic energy (this was before 
the QPOs were discovered.) 

The exceptional energetics of the SGR 1806 — 20 
hyperflare Ij, fl4| . the close proximity of the source 
pi l28l. I29I [30j | and the availability ofprecisely measured 
QPO frequencies and bandwidths [a 0, @] made SGR 
1806 — 20 attractive for study as a possible GW emitter. 

In this paper we make use of the LIGO Hanford (WA) 
4 km detector (HI), the only LIGO detector collecting 
low noise data at the time of the flare, to search for or 
to place an upper bound on the GW emission associated 
with the observed QPO phenomena of SGR 1806 — 20. 
At the time of the event the GEO600 detector was also 
collecting data. However, due to its significantly lower 
sensitivity at the frequencies of interest, it was not used 
in this analysis. 

As will be shown, the 92.5 Hz QPO upper bounds can 
be cast into a characteristic GW energy release in the ~ 
8xl0 46 - 3xl0 47 erg (~ 4xl0~ 8 - 2xlO- 7 M c 2 ) range. 
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Observation Frequency FWHM Period Satellite References 

[Hz] [Hz] [s] 

a 17.9 ± 0.1 1.9 ± 0.2 60-230 RHESSI [7] 

b 25.7 ± 0.1 3.0 ± 0.2 60-230 RHESSI [7] 

c 29.0 ± 0.4 4.1 ± 0.5 190-260 RXTE [8] 

d 92.5 ± 0.2 1.7+o 4 170-220 RXTE [6] 

e " " 150-260 " [8]" 

f 92.7 ± 0.1 2.3 ± 0.2 150-260 RHESSI [7] 

g 92.9 ± 0.2 2.4 ± 0.3 190-260 RXTE [8] 

h 150.3 ± 1.6 17 ± 5 10-350 RXTE [8] 



i 



626.46 ± 0.02 0.8 ± 0.1 50-200 RHESSI [7] 



1 625.5 ± 0.2 1.8 ± 0.4 190-260 RXTE 

m 1837 ± 0.8 4.7 ± 1.2 230-245 RXTE 



a Ref. Q makes an adjustment to the observation period of Ref. 

a " 

TABLE I: Summary of the most significant QPOs observed in the pulsating tail of SGR 1806 — 20 during the 27 December 
2004 hyperflare (from Ref. [j|). The period of observation for the QPO transient is measured with respect to the flare peak, 
the frequencies are given from the Lorenzian fits of the data and the width corresponds to the Full-Width-at-Half-Maximum 
(FWHM) of the given QPO band. 



This energy approaches the total energy emitted in the 
electromagnetic spectrum and offers the opportunity to 
explore the energy reservoir of the source. In the event 
of a similar Galactic hyperflare coinciding with LIGO's 
fifth science run (S5), the energy sensitivity involved at 
- 100 Hz would probe the - 2 x 10 45 erg (~ 10" 9 M c 2 ) 
regime. 



cant components at ~ 150 and ~ 1840 Hz as well. Table 
U is taken from Ref. [8J] and summarizes the properties 
of the most significant QPOs detected in the X-ray light 
curve tail of the SGR 1806 - 20 giant flare. 



III. THE LIGO DETECTORS 



II. SATELLITE OBSERVATIONS 

SGR 1806 — 20 is a Galactic X-ray star thought to be 
at a distance in the 6 to 15 kpc range [J, [H, [2^, [30| . 
The total (isotropic) electromagnetic flare energy for the 
27 December 2004 record flare was measured to be ~ 
10 46 ergs [1, [HI assuming a distance of 10 kpc. 

QPOs in the pulsating tail of the SGR 1806-20 hyper- 
flare were first observed by Israel et ai.[(J using RXTE, 
and revealed oscillations centered at ~ 18, ~ 30 and 
- 92.5 Hz. Using RHESSI, Watts and Strohmayer 
confirmed the QPO observations of Israel et al. reveal- 
ing additional frequencies at ~ 26 Hz and ~ 626.5 Hz 
associated with a different rotational phase. Closer in- 
spection of the RXTE data by Strohmayer and Watts Q 
revealed a richer presence of QPOs, identifying signifi- 



Thc Laser Interferometer Gravitational Wave Observa- 
tory (LIGO) [12] consists of three detectors, two located 
at Hanford, WA (referred to as HI and H2) and a third 
located in Livingston, LA (referred to as LI). Each of 
the detectors consists of a long-baseline interferometer 
in a Michelson configuration with Fabry-Perot arms (see 
Ref. [3l| for details). The passage of a GW induces a 
differential arm length change AL which is converted to 
a photocurrent by a photosensitive element monitoring 
the interference pattern of the detector. This electrical 
signal is then amplified, filtered and digitized at a rate of 
16384 Hz to produce a time series which we refer to as 
the GW channel. 

To calibrate the GW channel in physical units, the in- 
terferometer response function is frequently measured by 
generating known differential arm length changes. The 
uninterrupted monitoring of the response function is en- 
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sured with the addition of continuous sinusoidal excita- 
tions referred to as calibration lines. 

The interferometer sensitivity to AL enables to mea- 
sure a strain h defined as 



where L denote the mean of the two arm lengths. The 
target frequency range of interest is the audio band with 
frequencies in the 50 Hz to 7 kHz range. 

LIGO has dedicated science runs when good and re- 
liable coincidence data is available, alternating with pe- 
riods of commissioning to improve the sensitivity of the 
instrument. In order to cover times when an astrophys- 
ically notable event might occur, such as the 27 Decem- 
ber 2004 event of this analysis, data from times when 
commissioning activities do not disable the machine is 
archived by a program referred to as Astrowatch [32| . 
Due to the nature of the time period, the detector's con- 
figuration was continuously evolving and was not as well 
characterized as the dedicated science runs. On the other 
hand, there was a deliberate attempt to place the inter- 
ferometers in a high-sensitivity configuration compatible 
with the commissioning modifications of the epoch. 

At the time of this event two of the LIGO detectors 
were undergoing commissioning in preparation for the 
fourth science run (S4). Only data from HI is available 
for the analysis of this event. 

Figure Q] plots the best strain-equivalent noise spectra 
of HI during the S4 and S5 data-taking periods (light 
gray curves). The average noise spectra at the time of 
the flare is shown by the dark gray curve and the dashed 
line describes the design sensitivity. 



IV. DATA ANALYSIS 

This analysis relies on an excess power sear ch [43l| , vari- 
ants of which are described in Ref. [HI, [U, l35j|. In this 
analysis we compare time-frequency slices at the time 
of the observations with neighboring ones. The algo- 
rithm used analyzes a single data stream at multiple fre- 
quency bands and can easily be expanded to handle coin- 
cident data streams from multiple detectors. The trigger 
provided for the analysis corresponds to the flare's X- 
ray peak as provided by the GRB Coordinate Network 
(GCN) reports 2920 0| and 2936 [H[ at time corre- 
sponding to 21:30:26.65 UTC of 2004-12-27. 

In the absence of reliable theoretical models of GW 
emission from magnetars, we keep the GW search as 
broad and sensitive as possible. The search follows the 
QPO signatures observed in the electromagnetic spec- 
trum both in frequency and time interval. In particular, 
we measure the power (in terms of detector strain) for 
the intervals at the observed QPO frequencies (as shown 
in Tab. [TJ for a given bandwidth (typically 10 Hz) and we 
compare it to the power measured in adjacent frequency 



bands not related to the QPO. The excess power is then 
calculated for each time-frequency volume of interest. 

Although QPOs are not observed in X-rays until some 
time after the flare, the magnetar model suggests that 
the seismic modes would be excited at the time of the 
flare itself. For this reason, we also search for GW emis- 
sion associated with the proposed seismic modes from 
the received trigger time of the event. In addition, we 
chose to examine arbitrary selected frequency bands, re- 
ferred to as control bands, whose center frequency is set 
to twice the QPO frequency and processed identically to 
the QPO bands. This allowed us to cover a wider range 
of the detector's sensitivity while allowing the reader the 
flexibility to estimate the sensitivity to low significance 
QPOs not addressed here (see Ref. @) as well as future 
observations/exotic models of GW emission yet to come. 

Another aspect of the satellite observations is the 
quasi-periodic nature of the emitted electromagnetic 
waveform with a possible slow drift in frequency. Since 
there is no knowledge of the GW waveforms that would 
be associated with this type of event, we tune our search 
algorithm to be most sensitive to long quasi-periodic 
waveforms with fairly narrow bandwidths while short 
bursts are strongly discriminated against. The waveform 
set used in testing the sensitivity of the algorithm by 
adding simulated data in the analysis software is chosen 
in line with this argument. 



A. Pipeline 

A block diagram of the analysis pipeline is shown in 
Fig. [3 where the Gamma-ray bursts Coordinates Net- 
work (GCN) reports provide the trigger for the analy- 
sis. The on and o/f-source data regions are then selected 
where the former corresponds to the QPO observation 
periods, as shown in Tab. HI The off-source data region 
begins at the end of the six minute long QPO tail (set to 
400 s after the flare peak) lasting to ten minutes prior to 
the end of the stable HI lock stretch for a total of ~ 2 h 
of data. 

The on-source region consists of a single segment. 
This segment either starts at the moment of the flare 
(istart = to) or at the beginning of the QPO observation 
(^start = ^qpo) and lasts until the end of the observation 
(tend)' The off-source region consists of numerous non- 
overlapping segments, each of duration At = i cn d — istart ■ 

To provide an estimate of the search sensitivity, an ar- 
bitrary simulated gravitational waveform can be added 
(or injected) to each off-source data segment. All of the 
segments ( on- or ojf-source) are processed identically. In 
the procedure described by the conditioning block, the 
data is band-pass filtered to select the three frequency 
bands of interest: the QPO band as shown in Tab. |T] 
and the two adjacent frequency bands. Using the inter- 
ferometer response function at the time of the event, the 
data is calibrated into units of strain and a data-quality 
procedure, as described below, is applied to the data set. 
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FIG. 1: The strain-equivalent sensitivity of the HI detector at the time of the hyperfiare, the fourth and fifth science runs (S4, 
S5), and its design sensitivity. 
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FIG. 2: A block diagram of the analysis sketching the signal 
flow. 



After the conditioning procedure is complete, the data 
stream is pushed through the search algorithm, which 
computes the power in each segment for the three fre- 
quency bands of interest and then the excess power in 
the segment. Finally, on- and off-source excesses are com- 
pared and in the case of no significant on-source signals, 
the standard Feldman-Cousins [36| statistical approach is 
used to place an upper limit based on the loudest signal. 



The data processing can be validated against analyt- 
ical expectations by replacing the off-source region with 
simulated data. 



B. Data conditioning 

The conditioning procedure consists of zero-phase fil- 
tering of the data with three different band-pass Butter- 
worth filters. The first band-pass filters the data around 
the QPO frequency of interest with a predefined band- 
width. This bandwidth depends on the observed QPO 
width (see Tab. QJ and on the fact that the QPOs have 
been observed to evolve in frequency. For the QPOs ad- 
dressed here, the bandwidth is set to 10 Hz (well above 
the measured FWHM shown in Tab. [I] with the excep- 
tion of the 150.3 Hz oscillation where the bandwidth was 
set to the measured FWHM, 17 Hz. 

The bandwidth for the control bands is also set to 
10 Hz which is still above twice the measured FWHM. 
An exception to this is the 150.3 Hz second harmonic 
which is within one Hz away from the fifth harmonic of 
the 60 Hz power line. The bandwidth in this case is set 
to twice the measured FWHM (2 * 17 Hz = 34 Hz) but 
a 4 Hz wide notch at 300 Hz is included to suppress the 
significant sensitivity degradation provided by the line. 
For this reason, the effective bandwidth is 30 Hz. 

The data is also filtered to select the two adjacent fre- 
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quency bands with identical bandwidths of the chosen 
QPO band. Using the adjacent frequency bands allows 
us to discriminate against common non-stationary broad- 
band noise, thereby increasing the search sensitivity, as 
will be described in Sec. IIVC1 

A gap between frequency bands was introduced for 
some of the QPO frequencies in order to minimize the 
power contribution of known instrumental lines. Fur- 
thermore, 60 Hz harmonics which landed in the bands 
of interest were strongly suppressed using narrow notch 
filters. 

The three data streams are calibrated in units of strain 
using a transfer function which describes the interferom- 
eter response to a differential arm length change. 

The conditioning procedure ends with the identifi- 
cation of periods of significant sensitivity degradation. 
These periods are selected by monitoring the power in 
each of the three frequency bands in data segment dura- 
tions, or tiles, 125 ms and 1 s long. If the power is above a 
set threshold in any of the three bands, the tile in ques- 
tion identifies a period of noise increase. This abrupt 
power change in a second-long time frame (or less) does 
not correspond to a GW candidate lasting tens to hun- 
dreds of seconds long. For this reason, the full data set 
contained in the identified tile is disregarded and short- 
duration GW bursts, not among the targeted signals, 
would be excluded by this analysis. 

To set a particular threshold we first determined the 
variance of the resulting power distribution which was 
calculated by removing outliers iteratively. As will be 
described in Sec. |V] we used 2a, 3a and 4a cuts and 
we injected different waveform families to optimize the 
search sensitivity. 

C. The search algorithm 

The algorithm at the root of the search consists of tak- 
ing the difference in power between a band centered at 
a frequency / qpo and the average of the two frequency 
bands adjacent to the QPO frequency band, also of band- 
width A/, typically centered at f± = / qpo ± A/. 

After band-pass filtering, we are left with three chan- 
nels for each QPO: c qpo (t), c + (t), and c_(i). The power 
for the QPO interval is for each of these channels: 

^qpo,±= /" ld (c qpo ,±) 2 ^ (2) 

•/Wrt 

where tiles that were vetoed are excluded from the inte- 
gral. The excess power is then defined as 

AV = V qpo - 7> avg (3) 

where P avg = (V++V-)/2 is the average of the adjacent 
bands. We refer to the resulting set of AV calculated 
over the off-source region as the background while the on- 
source region provides a single excess power measurement 
of duration At for the period from t s tart to tend- 



V. SENSITIVITY OF THE SEARCH 

In order to estimate the sensitivity of the search, differ- 
ent sets of more or less astrophysically-motivated wave- 
forms, or in some cases completely ad-hoc waveforms, are 
injected in the off-source region and the resulting excess 
power is computed. 

The strength of the injected strain (at the detector) 
hdct{t) is defined by its root-sum-square (rss) amplitude, 
or 

/ /-ti+At 

^rss-dot = J J |/ldct(i)| 2 dt (4) 

integrated over the interval At, as described in Sec. IIV A[ 
where t\ indicates the start of a segment in the back- 
ground region. The search sensitivity to a particular 
waveform, ^s^-det' ^ s defined as the injected amplitude 
'irss-dct such that 90% of the resulting AT is above the 
off-source median. This choice of definition provides 
a characteristic waveform strength which, on average, 
should not be far from a 90% upper bound. 

We injected various waveform families (namely Sine- 
Gaussians (SG), White Noise Bursts (WNB), Amplitude 
(AM) and Phase Modulated (PM) waveforms) in the off- 
source region to quantify the sensitivity of the search to 
these types of waveforms. Each waveform was added di- 
rectly to the raw data segments and the search sensitivity 
was explored as a function of the various parameters. As 
previously mentioned, we designed the algorithm to be 
sensitive to arbitrary waveforms with a preset small fre- 
quency range while discriminating against any type of 
short duration signals. 

The result of the sensitivity study for the case of the 
92.5 Hz QPO (observation d of Tab. P) is shown in Fig. [3] 
where the band center frequencies, bandwidths and signal 
durations were set to / qpo = 92.5 Hz, /_ = 82.5 Hz, 
/+ = 102.5 Hz, Af = 10 Hz and At = 50 s. 

SG waveforms are parameterized as follows 

h dct (t) = A sin (2nf c t + <j>) e -(*-*°) 2 /r 2 (5) 

where A is the waveform peak amplitude, f c is the wave- 
form central frequency, Q — V2iTTf c is the quality factor, 
r is the 1/e decay time, is an arbitrary phase and to 
indicates the waveform peak time. In the case of Q — ► oo 
the waveform approaches the form of a pure sinusoid. 
The top left panel of Fig. [3] plots the search sensitivity 
versus the quality factor Q of the injected SG waveform, 
indicating that the analysis is most sensitive to SG wave- 
forms with quality factors in the range Q € [~ 10 3 : oo]. 
The response is also shown as a function of a 2a and 
4a data quality cut on the off-source RMS distribution 
calculated for 125 ms long tiles. The more aggressive 
2a cut yields significantly better results and was chosen 
for the 92.5 Hz QPO analysis. This band in particular is 
significantly more problematic than the others exhibiting 
a high-degree of non-stationarity as well as a relatively 
high glitch rate. 
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FIG. 3: Search sensitivity to different waveform families and for different data quality cuts. The cuts are relative to the off- 
source RMS distribution calculated in segments 125 ms long and for 2er cuts (dark gray crosses) and 4<r cuts (light gray crosses). 
Top left: SG waveforms injections as a function of quality factor Q varied from Q = 600 to Q = 10 6 . Dashed line represents the 
average sensitivity (5.1 x 10~ 22 strain Hz~ 1//2 ) for injections with Q > 5 x 10 3 (where the sensitivity is essentially flat) and a 
2cr cut. Top right: 40 s long WNBs waveform injections as a function of burst bandwidth ranging from 1 Hz to 11 Hz. Within 
the parameter space explored the sensitivity is essentially constant. Bottom left and right: PM and AM waveform injections 
as a function of modulation depth for a modulation frequency of 100 mHz. 



The decline in sensitivity as the Q decreases originates 
from the data quality procedure. As parameter Q takes 
smaller values, the waveform energy begins to concen- 
trate in shorter time scales and the conditioning pro- 
cedure identifies and removes intervals of the injection 
which are above threshold. In the 2a case, the sensitiv- 
ity is relatively flat for Q > 5 x 10 3 and the average value 
is ft-rss-det = 5-1 x 10~ 22 strain Hz -1 / 2 also shown in the 
plot by the dashed line. The corresponding waveform 
duration 5t, defined as the interval for which the wave- 
form amplitude is above A/e, is St > \[2 Q/irf c ~ 24 s, 
appropriate for the targeted search as shown in Tab. [U 

The top right panel of Fig. [3] plots the sensitivity to a 
large population of 40 s long WNBs injections of band- 
widths ranging from 1 Hz to 11 Hz. The waveform is gen- 
erated by band-passing white noise through a 2 nd order 
Butterworth filter with bandwidth defined at the -3dB 
cutoff point and burst duration set by a Tukey window. 
As shown in the SG case, the most aggressive 2er cut 
outperforms the 4cr and no significant departure in sen- 
sitivity is seen for bandwidths up to 10 Hz. It is worth 
noting that WNBs would correspond to incoherent mo- 
tion of the source and may not be physical. However the 



purpose of this study is to quantify the robustness of the 
search to a variety of waveforms. 

The bottom two panels of Fig. [3] plot the sensitivity to 
PM and AM waveforms versus modulation depth, where 
the modulation frequency is set to / mo d = 100 mHz for 
both cases. These waveforms are used to investigate QPO 
amplitude and frequency evolutions. For the PM case, 
the waveform is described as 

hdctif) = A cos (27r/ c t + k mod x(t) + (jy) (6) 

where A is the waveform amplitude, / c is the carrier fre- 
quency, 4> is an arbitrary phase, A; mo d is a modulation 
depth constant and x(t) is the modulation signal 

x[t) = sin (2n f mod t) (7) 

It can be shown that the instantaneous frequency / is 

f(t) =f c + A/ mod cos(27r/ mod t) (8) 

where A/ mod = fc mod / mod . From Fig. [3] the PM sen- 
sitivity is essentially constant within modulation depths 
in the range A/ mod € [1 : 5] Hz. 
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The AM injection is parameterized as 



h det (t) = A(t) cos (2nf c t) 



where 



A{t) = A 



sin(27r/ mo di) - "'mod 



1 



(9) 



(10) 



with waveform constant amplitude Ao, fcmod modulation 
constant, and f c carrier frequency. The search sensitivity 
to this waveform family can be expressed in terms of the 
modulation depth R defined as 



R = 1 



1 - k 



mod 



l 



v mod 



i 



'•'mod 



(ii) 



The bottom right panel of Fig. [3] plots the sensitivity 
of this waveform as a function of R. As fc mo d — * oo, the 
modulation depth parameter R — > 0, no modulation is 
applied and the waveform is a sinusoid of constant am- 
plitude. As /c mo d — > 1, the modulation depth is maximal 
(R = 1) and the amplitude A(t) is also sinusoidal in na- 
ture. From Fig. [3] the AM sensitivity is essentially con- 
stant within modulation depths in the range R € [0 : 1]. 
The average response to SG, as shown in the top left 
panel of Fig. [3j is also shown in the other three panels 
for comparison. 

The results shown in Fig. [3] indicate that the search 
sensitivity is approximately the same for all the wave- 
forms considered. 

It is also possible to estimate the theoretical search sen- 
sitivity to a sinusoidal injection. Assuming white gaus- 
sian stationary noise for the detector output, we can de- 
rive (see Ref. |35[ ) the following expression for the search 
sensitivity 



-.thco 



s— det 



1.25 Sl'\f) (Af At) 



1/4 



(12) 



1 /2 

where SJ (/) is the strain-equivalent amplitude spectral 
density of the detector noise at frequency /, in units of 
strain Hz^ 1 / 2 , and Af and At are the bandwidth and 
duration of the segment in question, in units of Hz and s. 
The order-of-unity factor (1.25) stems from the 90% sen- 
sitivity definition as previously discussed and from taking 
the difference in power between bands. 

Referring to Fig. [TJ the strain sensitivity at / = 
92.5 Hz is Sl /2 (f) ~ 9 x 10~ 23 strain Hz" 1/2 . Using 
Af = 10 Hz and At = 50 s, the expected sensitivity is 



- theo 
rss— det 



5.3 x 10" 22 strain Hz 



-1/2 



(13) 



in good agreement with the average response of 



L sens 
"rss-dct 



5.1 x 10~ 22 strain Hz 1/2 shown in Fig. [31 



VI. RESULTS 

Inspection of the on-source data segments revealed no 
significant departure from the off-source distribution and 



we cast the results of this analysis in terms of upper 
bounds on GW signals. These limits are found to be 
well below the maximum allowed upper bounds in the 
non-detection regime, which we refer to as non-detection 
threshold, assuming a continuous observation of SGR 
1806 — 20 and requiring an accidental rate of one event 
in one-hundred years (see Tab. HI)) . 

We used the unified approach of Feldman-Cousins (36| , 
which provides upper confidence limits for null results, 
two-sided confidence intervals for non-null results and 
treats confidence limits with constraints on a physical re- 
gion. In view of the fact that at the time of the hyperflare 
event only one of the three LIGO detectors was collecting 
data and that the full detector diagnostic capability was 
not fully exploited, the lower bounds on the confidence 
intervals was set to zero (i.e. no detection claim based 
purely on the statistical analysis was allowed). 

Table X of Ref. [3(| was used to place the upper limits 
of this search. The excess power distribution for the off- 
source region of each QPO transient was parameterized 
with a Gaussian Probability Density function (PDF), and 
the mean jj,, standard deviation a and their relative errors 
is estimated. The on-source excess power measure and 
the lookup table were then used to set 90% confidence 
intervals. 

Table [H] presents the results of this search, for both 
the control and QPO frequencies, in terms of 90% upper 
bounds on the GW waveform strength, h: 



90% 
rss— det ' 



sured at the time of the observation. The first column 
of the table indicates the observation we address, with 
reference to the original measurements shown in Tab. HI 
The second, third, fourth and fifth columns indicate the 
center frequency, bandwidth, period, and duration used 
in the search. The sixth column, labeled as non-detection 
threshold, lists the maximum upper bound allowed in the 
non-detection regime. A data quality flag was used for 
the 92.5 Hz QPO observation only, with a power thresh- 
old set at the 2a level relative to tiles 125 ms long. 

The last column, labeled h^° s ^ dct , presents the results 
where the contributions due to the different uncertainties 
are shown separately. The first of these, the first number 
in superscript, shows the 90% upper bound arising from 
the statistical uncertainties in the off-source estimation. 
These uncertainties are generated using a Monte-Carlo 
simulation: a set of means fx and standard deviations 
a are extracted from Gaussian distributed populations 
of standard deviation ap, and <j& corresponding to the 
fit parameter uncertainties. For each (/t, a) combination 
and the same on-source excess power measure we used 
the lookup table in Ref. [36j to generate 90% confidence 
intervals for the quoted upper limit. 

The second uncertainty quoted is statistical and arises 
from errors in the detector response function to GW ra- 
diation via the calibration procedure. We placed a con- 
servative estimate of the calibration accuracy to a one 
standard deviation of 20%. The third uncertainty is a 
systematic error of 6% also arising from the calibration 
procedure. 
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Observation Frequency Bandwidth Interval Duration 
[Hz] [Hz] [s] [s] 


Thresholdnon-det 
[lO^strain Hz" 1/2 ] 


v90% 
rss — det 

[10- 22 strain Hz~ 1/2 ] 


e,f 


92.5 


10 


150-260 


110 


lo.U 


2.75 


+0.47 +0.70 +0.16 +0.77 
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4.53 


g 
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2.90 


+0.43 +0.74 +0.17 +0.75 





4.67 


(1 
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9.50 


control freq. 


185.0 


8 


150-260 


110 


ion 


9.48 


_i_n j_9 zi^ j_n ^7 _i_n 97 

-f-U.OJ- ^^.*±0 T^.u' T'J.Z/ 




12.8 








190-260 


70 


17.6 


8.17 


+0.40 +2.09 +0.49 +0.17 


= 


11.0 








170-220 


50 


10. D 


8.03 


i o on _lo nfi j_n aq. _i_n oa 
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24.1 


11.4 


+ 1.06 +2.91 +0.68 +0.00 




15.1 


h 
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17 


0-350 
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30.2 


12.4 


+ 1.78 +3.16 +0.74 +0.00 


= 


16.7 


control freq. 


300.6 


30 


0-350 
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26.4 


+4.46 +6.75 +1.58 +0.00 
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+ 1.76 +6.56 +1.54 +0.00 
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1 
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19.4 


-Li 9^ j-A Q7 j-1 17 _i_n nn 
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0-260 


260 


60.1 


28.2 


+2.70 +7.22 +1.69 +0.00 




37.6 


control freq. 


1253.0 


10 


50-200 


150 


114 


49.4 


+4.10 +12.64 +2.96 +0.00 




65.6 








190-260 


70 


89.0 


30.6 


+2.69 +7.84 +1.84 +0.00 




40.7 








0-260 


260 


107 


53.5 


+4.50 +13.71 +3.21 +0.00 




71.2 


m 


1837.0 


10 


230-245 


15 


94.7 


34.6 


+ 1.26 +8.86 +2.08 +0.00 




45.6 








0-245 


245 


192 


54.9 


+ 11.72 +14.05 +3.29 +0.00 




76.5 



TABLE II: List of frequencies and observation times used in this analysis with the corresponding results. The first column 
describes the addressed QPO observation, labeled by letters as they appear in Tab. U A wider range of the detector's sensitivity 
can be explored using the frequency bands here labeled as control frequencies (see text). The second, third, fourth and fifth 
columns indicate the center frequency, bandwidth, interval, and duration used in the search. The sixth column provides the 
non-detection threshold. The last column presents the results where the contributions due to the different uncertainties are 
shown separately. The first two numbers in superscript represent the statistical uncertainty in the off-source estimation and 
calibration procedure respectively. The third one shows the contribution of a systematic uncertainty of 6% due to the calibration 
procedure. The last uncertainty is a systematic arising from the off-source data modeling which depends on the presence of 
outliers (see text for details). To produce the upper bound h*^ dct statistical contributions are added in quadrature while the 
systematic contributions are added linearly. 



The occasional presence of tails in the off-source seg- 
ments, consisting typically of a few large excess power 
measurements in the off-source data of each QPO intro- 
duces a bias in the upper bounds which is presented as 
a source of systematic uncertainty (represented by the 
fourth number in superscript). This bias is quantified by 
including and excluding the off-source distribution ±3<r 
outliers from the fitting procedure and the difference in 
the upper bounds, SHg?^ = K^ dct -h^° d f t is shown 
in the column in question. 

In order to fold in the different uncertainties we sum in 
quadrature the statistical uncertainties shown (originat- 
ing from the off-source estimation and the calibration) 
and we increase the bound by the two systematic errors. 



VII. ASTROPHYSICAL INTERPRETATION 

In this section we provide a characteristic GW en- 
ergy -Eqw associated with the measured upper bounds 
^rss-det' snown m Tab. [Ill cast in terms of a simple 
source model. In this model we assume that the emission 
is isotropic, that the plus and cross polarization states are 
uncorrelated but have equal power. 

Under these assumptions (equal uncorrelated power ra- 
diated in the plus and cross polarizations) the strain in 
the detector can be related to the GW flux incident on 
the Earth via 

h 2 Tas - de t = \(Fl+F*)hi s (14) 
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where 



hL = 



K(t) 



h 2 x (t)]dt 



(15) 



and F+ and F x are antenna response functions that de- 
pend on (i) the right-ascension and declination of the 
source, (ii) the time of the flare, (hi) the location and 
orientation of the detector, and (iv) a polarization an- 
gle defining the plus and cross polarizations. The depen- 
dence on this polarization angle vanishes in the combina- 
tion F++F%, which is a quantity ranging from to 1; the 
Hanford detector's antenna response to SGR 1806 — 20 
at the time of the hyperflare was 



Fl + Fl = 0.174 



(16) 



This shows that the source was not particularly well sit- 
uated in the detectors antenna pattern. Under our as- 
sumption of isotropic emission, the energy released by 
the source is related to the gravitational wave flux at the 
Earth by 



^2 3 2 f 2 
fi GW — ~ ''">■ 



G 



(17) 



In terms of the upper limits presented, the equivalent 
bound on the gravitational wave emission corresponding 
to a particular QPO is 



E^ 90% = 4.29 x 10- 8 M Q c 2 x 



(18) 



2 / f x 2 



u90% 

rss — det 



lOkpc/ V92-5 Hz / I 4.53 x 10 - 22 strain Hz" 1/2 



(here the values of the best QPO strain bound are used). 
It is worth noting that the best energy upper bound is 
comparable to the energy emitted in the electromagnetic 
spectrum (see for example Ref. Q). 



VIII. CONCLUSION 

Quasi-Periodic Oscillations have been observed in the 
pulsating X-ray tail of the SGR 1806 — 20 hyperflare of 
27 December 2004 by the RXTE and RHESSI satellites. 
The present consensus interprets the event as a dramatic 
re-configuration of the star's crust and/or magnetic field. 
In turn, this starquake could plausibly excite the star's 
global seismic modes and the observed QPOs could po- 
tentially be driven by the seismic modes. The energetics 
of the event, the close proximity of the source, and the 
availability of observed QPO frequencies and bandwidths 
provided a unique opportunity to measure GWs associ- 
ated with this phenomenon. 

Upper limits in the gamma and high-energy neutrino 
flux were recently measured by the AMANDA-II detec- 
tor [37|]. However the only other published GW search 
associated with the SGR 1806 — 20 hyperflare used the 
AURIGA bar detector [38| to place upper limits on the 



GW waveform strength emitted for frequencies around 
~ 900 Hz. At the time of the event, Hi's strain noise 
equivalent in the ~ 900 Hz region is a factor ~ 5 lower 
than AURIGA'S. 

The AURIGA search targeted different physics, there- 
fore the comparison to our results is not possible. Expo- 
nentially decaying sinusoids of decay time 100 ms were 
searched for by measuring the power in time and fre- 
quency slices of At = 201.5 ms and A/ = 5 Hz respec- 
tively in the 855 Hz to 945 Hz range. A set of 95% 
upper bounds on the waveform strength were placed in 
the \ 9 s f° dct = 1.4 x 10- 21 strain Hz" 1 / 2 to h»™ det = 
3.5 x 10~ 21 strain Hz~ 1/2 range. 

At the time of the event one of the three LIGO de- 
tectors was in operation under the Astrowatch program. 
Under this program, data is collected at times of com- 
missioning when the interferometers are not undergoing 
adjustments. Only ~ 2 h of data was available for this 
analysis. 

An algorithm was designed to measure the excess 
power deposited in the machine at the time of the 
event. This algorithm exploits power measures in mul- 
tiple bands to reject common mode noise sources, such 
as broadband noise. Power measures in time scales less 
than Is are also monitored to reject fast signatures in- 
consistent with the scope of this analysis. 

The design was driven by the desire to repeat this mea- 
surement for future flares with the ability to use multiple 
data streams from multiple detectors, focusing on mod- 
ularity, flexibility, and simplicity. 

Signals were software-injected into the raw data stream 
to study the analysis sensitivity to a variety of waveform 
families and parameters. A large astrophysical motivated 
parameter space was explored under which the search 
sensitivity is essentially constant. 

At the time of the event, the strain-equivalent ampli- 
tude spectral density of the detector output was a fac- 
tor of a few away from the one corresponding to the 
fourth science run. Under this condition, the best up- 
per limit that we place corresponds to the 92.5 Hz QPO 
observed 150 s to 260 s seconds after the flare. In terms 
of waveform strength, we place a 90% upper bound of 
h 90% = 4 53 x 1Q -22 strain Hz -i/2 on the Gw wave _ 

rss— QGt 

form strength in the detectable polarization state reach- 
ing our Hanford (WA) detector, which, in terms of a 
simple source model, provides a characteristic energy 

This 



«f U7 ° = 7.67 x 10 46 erg (4.29 x 10" 8 M Q c 2 ) 



r iso,90% 
/GW 

is the best upper limit published on the GW waveform 
strength on this type of source and represents the first 
multiple-frequency asteroseismology measurement using 
a GW detector. It is also worth noting that this en- 
ergy estimate is of the same order as the isotropic energy 
estimate measured electromagnetically, providing the op- 
portunity to probe the energy reservoir of the source. 

The limits presented here represent GW strength ob- 
tained by the LIGO detectors in late 2004. At the 
time of this writing, LIGO is undergoing a data-taking 
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period, referred to as the fifth science run S5, where 
all three interferometers have reached design sensitiv- 
ity, [3§|. The improvement at 150 Hz corresponds to 
a decrease in strain-equivalent noise of > 3 in terms of 
GW energetics. This estimate excludes the sensitivity 
increase that can be achieved by cross-correlating data 
streams from the multiple LIGO detectors. A follow-up 
of this analysis will certainly examine the various SGR 
1806 - 20/SGR 1900 + 14 outbursts, which occurred in 
the 2005 - 2006 period, exploring GW energetics which 
probe the ~ 2 x 10 45 erg (~ 10~ 9 M Q c 2 ) regime. 

At the end of the S5 data-taking period, the initial 
LIGO detectors will be upgraded to an enhanced state 
[40l | which we refer to as Enhanced LIGO. The foreseen 
improvement will be a factor of ~ 2 in strain-equivalent 
noise for frequencies above 100 Hz. The future GEO- 
HF [4l[ detector will provide a significant high-frequency 
improvement in sensitivity providing an opportunity to 
study future high-frequency QPOs. 

Advanced LIGO [46[ will provide an increase in strain- 
equivalent sensitivity of ~ 10 with respect to the initial 
LIGO detectors while opening up the low (10 — 50 Hz) 
frequency range. This offers a particularly interesting 
opportunity because a lower frequency search would be 
feasible. For hyperflare events occurring at the time of 
its operation, the observable GW energetics at 100 Hz 
would lie in the - 2 x 10 43 erg (~ 10~ n M Q c 2 ) regime. 
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